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Procede et dispositif d'^galisation non liniaire d'un signal multi-porteuse dans une liaison satellite. 
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fe?) Le proc6d6 consisle d effectuer (1) une post-distorslon 
dusignal provenant du satellite, d decoder (2) les points de 
la constellation du signal multi-porteuse regu et d corriger 
4) chaque point de la constellation par Elimination des 
^chos fixes pour estimer les points de la constellation qui 
sont les plus proches des points corripEs. 

Les applications vont k la realisation d'^gat'iseurs pour 
r^cepteurs de t^lEvision. 
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1 

Proc6d6 et dispositif d'^galisatton non linSaire 
d'un signal multi-porteuse dans une liaison satellite 

La pr6sente invention concerne un proc6d6 et un dispositif d'6gali- 
5 sation non Iin6aire d'un signal multi-porteuse d'une liaison satellite. 

Dans le cadre des transmissions de donn^es num6riques, diff6- 
rentes modulations peuvent etre envisag6es pour transmettre des infor- 
mations, selon le type de support de transmission utilis6. Dans le 
domaine hertzien les choix semblent actuellement plutot se porter vers 

10 une utilisation d'une modulation de type multi-porteuse. L'avantage 
obtenu est une uniformity de la modulation qui permet de transmettre le 
signal multi-porteuse sur des liaisons satellites. Cependant les non-lin6ari- 
t6s introduites par les tubes S ondes progressives qui 6quipent T^lectro- 
nique de retransmission des signaux des satellites imposent de pr6voir 

15 dans les stations terrestres de reception des dispositifs d'6galisation re- 
produisant la fonction de transfert inverse des tubes d ondes progres- 
sives. Malheureusement le caractdre hautement non-lln6aire de cette 
fonction de transfert rend tr6s difficile la reconstitution de la fonction de 
transfert inverse. En effet comme d^crit dans Particle de A.A. SALEH 

20 ayant pour titre "Frequency independent and frequency dependent 
nonlinear models of TWT amplifiers" publi6 dans la revue IEEE Trans, on 
Comm. vol. COM 29 Novembre 1981 pages 1715 ^ 1720, la fonction de 
transfert de type amplitude-amplitude et amplitude-phase d'un tube d 
ondes progressives peut etre assimil6e ^ une fonction Iin6aire sur une 

25 plage d6termin6e mais pour des amplitudes plus importantes il se produit 
un §crasement de Tamplitude et un d§phasage notable entre sa sortie et 
son entree. Le signal obtenu au niveau des r6cepteurs est par conse- 
quent fortement alt6r6 par ces distorsions. Le ph6nom§ne est encore 
plus aggrav6 avec un type de modulation multi-porteuses du type OFDM 

30 par exemple qui est I'abrSviation anglo-saxonne de "Orthogonal 
Frequency Division and Multiplexing" dont le spectre est semblable h 
celui d'un bruit blanc et poss^de comme on le sait une forte dynamique 
dans le domaine temporel. En fait, la plupart des modems disponibles 
pour les liaisons satellites utilisent des modulations monoporteuses et 
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principalement la modulation d quatre 6tats de phase QPSK. Les 6gali- 
seurs utilises par ce type de modulation sont d filtre transverse d decision 
directe ou bien d decision incorpor6e dans une boucle. lis ^quipent par 
example les modems connus sous les abr6viations 70 Mbits ANT, EF 
Data SDM 70, EF Data SDM 450 etc... Mais comme il s'agit de modula- 
tion mono-porteuse S faible dynamique dans le domaine temporel ces 
6galiseurs ne tiennent pas compte des non-lin6arit6s sp6cifiques des 
tubes d ondes progressives, bien qu'un §galiseur de ce type ait 6t6 pro- 
pose par M. J. PALICOT dans un article intitul6 "Egalisation de perturba- 
tions non-lin6aires" au Trelzl6me Colloque du GRETSI de Septembre 
1991. Mais I'efficacit6 de cette modulation reste limit^e du fait du faible 
nombre d'6tats de la constellation d'6tats caract6risant les modulations 
QPSK. 

Le but de I'invention est de pallier les inconv^nients pr6clt6s. 

A cet effet, I'invention a pour objet un proc6d6 d'6galisation non 
lin^aire d'un signal multi-porteuse dans une liaison satellite caract§ris6 en 
ce qu'il consiste : 

- d effectuer une post-distorsion du signal provenant du satellite 

- d decoder les points de la constellation du signal multi-porteuse 

regu 

et ^ corriger chaque point de la constellation par Elimination des 
6chos fixes pour estimer les points de la constellation qui sont les plus 
proches des points corrig^s. 

D'autres caract6ristiques et avantages de I'invention apparaltront 
ci-aprds S I'aide de la description qui suit faite en regard des dessins an- 
nexes qui repr^sentent : 

- la figure 1 les diffSrentes Stapes du procEdE d'6galisation selon 
I'invention mises sous la forme d'un organigramme ; 

- la figure 2 un premier exemple de mise en oeuvre d'un 6galiseur 

selon I'invention ; 

- la figure 3 un deuxi6me exemple de mise en oeuvre d'un 6gali- 

seur selon I'invention. 

Le proc6d6 selon I'invention qui est repr6sent6 par les Stapes 1 S 5 
sur la figure 1 consiste d effectuer selon une premiere 6tape une post- 
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distorsion des ^chantilions de signal regu provenant du satellite, d rScu- 
pgrer d T^tape 2 les points de ia constellation du signal muiti-porteuse 
repu, d corriger d V6tape 3 ciiaque point de la constellation par 
Elimination des 6chos fixes et d estimer d f §tape 4 les points de la 
5 constellation qui sont les plus proches des points corrigSs. 

La post-distorsion du signal regu qui est effectuSe d I'Etape 1 a 
lieu en prenant en consideration la r^ponse en amplitude du satellite et 
notamment celle de son tube a ondes progressives. La loi qui caract6rise 
cette r6ponse est d6crite dans Tarticle de A. A. Saleh pr6cit6. II s'agit 
10 d'une loi inversible. L'expression du module du signal r qui est appliqu6 d 
rentr6e du tube S ondes progressives en fonction de celui r' qui est 
obtenu par sa s ortie est de la forme 

r = -! — ^ . (1) 
r" 

Pour effectuer ia post-distorsion du signal repu le traitement 
15 consiste d effectuer une approximation de la caractSristique (r, 
4)) = f(r',<|>') oCi 0 et d^signent les phases respectives du signal appii- 
qu6 d TentrSe et du signal sortant du tube d ondes progressives par des 
gains complexes g^K 

En posant m'^' . exp(j. 8 g = m<^' exp(j.e<j') et en tenant 
20 compte du module et de la phase de g^iK la correction rcor* ^cor ®st 
obtenue en application des relations 

rcor-exp(j.0cor) = Q . r\ exp(j. <!)' ) 

= m«'.r'.exp(i(<|)' + e«>)) 

ce qui donne : 

rcor =m<i>.r' (2) 

et <f>cor = (t>' + e«> (3) 
25 Le gain gO) est appliqu§ au point si le module r" est sltu6 dans une 

plage D*^"* = [s'^*, s'^'^'"] oO S<j) d6signe le seuil de decision d'apparte- 
nance attach^ au gain g(i). 

Chacun des §chantilIons d'un m§me symbole de modulation multi- 
symboles OFDM est ainsi post-distordu- 
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L'^tape 2 qui consiste d rScup^rer les points de ia constellation de 
symboles regus s'effectue par un calcul d'une transform6e de Fourier 
rapide ou FFT sur un nombre d6ternnin6 d*6chantillons composant un 
mgnne symbole et qui permet en absence d*6chos, d condition que les 
5 gains g^J) aient 6t6 correctement calculus, de r6cup6rer les points de la 
constellation 6mis. Ces calculs peuvent par exemple etre effectu6s sur 
des symboles comportant 1024 6chantillons d I'aide d'une FFT d 1024 
points. 

En fait, bien que I'antenne de reception des r6cepteurs de satel- 
10 lites soit g6n6ralement conform6e pour etre trfes directive, il faut pour 
r6cup6rer les symboles 6liminer les ^chos proches qui pourraient nuire d 
la detection des symboles provenant du satellite. Cette Elimination peut 
§tre exScut^e d Taide par exemple d'un Egallseur fr^quentiel du type de 
celui qui est d§crit dans la demande de brevet PCT/FR 89 00546, et qui 
15 permet d'identifier la r^ponse fr^quentielle du canal de transmission entre 
rScepteur et satellite pour d^convoluer ensuite les paquets repus par 
cette r6ponse. 

L'identification a lieu par une Emission pEriodique d'un symbole 
particulier appel6 "paquet test" bien connu du r6cepteur. La transform6e 
20 de Fourier du paquet test repu, compar6e S celle du paquet test 6mis 
fournit alors raie par rale les coefficients de correction. 

Si {xn} d6signe les 6chantillons du signal 6mis, {x'n} d6signe les 
6chantillons du signal repu et {hn} d6signe la r6ponse impulsionnelle du 
canal alors x'n et Xn sont li6s par le produit de convolution 
25 {x'n} = {hn}*{xn} (4) 

Moyennant Tutilisation d"un intervalle de garde de dur6e suffi- 
samment grande pour absorber V6cho le plus grand, la relation (4) prend 
alors la forme d'un produit dans le domaine frfiquentiel 

X\ = Hk.Xk (5) 

30 ou {Xk} repr6sente les informations contenues dans les raies de 

symboles de la modulation multi-porteuse OFDM et {X'^} ces memes in- 
formations obtenues au r6cepteur, 

Le role des paquets test est d'estimer la r6ponse fr6quentielle du 
canal {H|c}. En d6signant par {X^tst} et {X'^^s*} les transform§es de 
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Fourier rapide des paquets tests S l'6metteur et au r^cepteur respective- 

ment, chacune des raies v6rlfie la relation 
v» tst 

Xk 

Les coefficients ainsi estim§s sont ensuite appllqu6s pour corriger 
chaque paquet de donn6es repu. 

Une fois que le signal a -616 distordu et corrigS par les coefficients 
de correction I'organe de decision estime le point de la constellation qui 
est le plus proche de chaque point corrig6. L'estimation a lieu t partir 
d'un crlt6re d'erreur d6fini d partir de I'erreur d'estimation entre chaque 
point corrig^ et le point estimS. 

En d6signant par Xn chaque 6chantlllon de signal 6mis par le satel- 
lite, par gtj) les gains partiels, par X jj* les produits des gains, par Vn les 
6chantillons de signal sur lesquels ont lieu le calcul de FFT, par les 
6chantillons de signal resultant du calcul de FFT, par les coefficients 
de correction, par Ck le produit de la correction, par §k l'estimation de 
chaque point de la constellation, par I'erreur d'estimation et par Ej le 
critfere d'erreur, le calcul d'estimation commence par associer d chaque 
6chantillon de signal Xn un produit de gain X ij* sulvant la relation : 
X!i' = X„.l^,j,(Xn) (7) 
I :{X) = 1 si X 6 D*j' 



avec 



•dHx) = Osi Xis D«' 



et 

chaque §chantillon Vn est ensuite obtenu en effectuant une somme des 
6chantillons X jj* pond6r6e par les gains partiels g<J) sulvant la relation : 

Yk = Sg*^'x<i! (8) 

Un calcul de FFT sur les 6chantillons Vn est ensuite r6alis6 suivant 
la relation : 

Yk = 71 sVo-Wn""^ (9) 
N n=0 

Les coefficients de correction Cj^ de chaque 6chantillon sont 
ensuite calcul6s suivant la relation ^ 
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C|c = Hk.Yk (10) 

L'erreur d'estimation ei^ est alors d6termin6e par la relation : 
ek=Ck-a|c (11) 

Le crit&re d'erreur est alors dSduit de ei^ par la relation : 
^ (12) 



k€l k 

dans laquelle d^signe un sous-ensemble des indices des raies 
Smises. 

Les gains partiels g^i) sont calcul6s en mettant en oeuvre un algo- 
rithme de type gradient. En assimiiant e,- d une fonction de cout telle 
que : 

er = F(G) (13) 

oO G d^signe un vecteur ayant pour composantes les parties r^elles et 
imaginaires de chacun des gains gO) et F est una fonction quadratique, le 
minimum de cette fonction de coOt conduit d une solution de la forme : 
Gt=ui = Gt=i - a.V Sr.G (14) 

En explicitant les d6riv§es partielles de par rapport aux parties 
r6el!es g et imaginaires g et en utilisant le symbole * pour designer 
des parties imaginaires conjugu^es. 
de r 
3^ 



kelk 



3(ck - a k) / X • ^ 



iT" 



(ck - a k) 



(15) 



= 2 Z Re 
keik 

= 2 2 Re 
keik L 



9(ck - ak) 



2 XRe 

2 Z 
kelkL 



dc k 
9^ 



©k 



(e 

HkN-1 



WN"^ek 



— X X I," . Wn . e k 



N 



rv> 0 



= 2 XRe[Hk FFT(k){x!i»}.ek*] 
kel k 

De maniSre analogue : 



(16) 



4 
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15 



= 2 ZRe 



(17) 



r 

kel k 



= 2 S Re 
kelk 



1% 



Z ^ (1) • . e k 



■] 



(18) 



= 2 ZRef^ Zi.x|i'.Wi'*.erl 
k€lk L J 

= 2 £Re[Hk i.FFT(k){xii»}.ek*] 
kelk 

En posant : 

A= 2 Z [ftk.FFT{k){xij>}.ek"] 
kelk 

Les dSrivies partielles calcul^es ci-dessus s'expriment alors : 
^.RelA, .19. 

La lol d'6volution des gain g<j) deviant ainsi : 



= Re(A) 



(20) 



= g|n-l)T<^* - a.{Re(A) + Re(jA)) 

= g(n-1)T"* - a. A* 
et fmalement : 



gin-DT = g(n-1)T''' " « 



2Hk.FFT(k){X<iJ}.ek* 
kelk 



(21) 



Ce qui pr6cdde fait apparaTtre des r^sultats interm6diaires qui sont 
les sorties de FFT calcul6es sur les 6chantillons avant la distorsion appor- 
t6e par les gains gO). Ces FFT s'introduisent de manlfere naturelle dans 
les calculs et modiflent la structure de I'^gaiiseur afin de faciliter son im- 
20 plantation. 
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Le proc6d6 d6crit pr6c6demment est realisable de la fapon repre- 
sentee par la structure simplifiee de r^galiseur de la figure 2. Cet egali- 
seur comporte un quantificateur 1 coupl6 S des amplificateurs 2-| d 2j de 
gain respectifs gO). Un circuit sommateur 3 realise la fonction (8) et ap- 
5 plique les 6chantillons yp calculus d un dispositlf de calcul FFT 4 par Tin- 
term6diaire d'un registre d d6calage 5. Les 6chantillons fournis par le 
dispositlf de calcul de FFT 4. sont appHquSs sur une preml6re entree 
d*operande d'un circuit multiplieur 6 qui regoit sur une deuxidme entree 
d'opSrande chacun des coefficients de correction Hi^ pour ^laborer les 

10 produits de correction C|^ calculus suivant la relation (10). Les valeurs 
des echantiilons C|^ sont ensuite quantifies pour un quantificateur 7 pour 
eiaborer les valeurs d 'estimation S^. 

Toutefols, la loi d 'evolution des coefficients g(^) suggere une 
seconde structure telle que celle representee d la figure 3, faisant appa- 

15 raTtre autant de FFT qu'il y a de niveaux de seuillage. Dans regaliseur 
repr6sente d la figure 3 les elements homologues ^ ceux de la figure 2 
sont repr6sent6s portent les memes references, A ia difference avec la 
figure 2 chacune des n sorties du circuit de quantification 1 est coupiee 
^ un circuit amplificateur 2^ ^ 2j par I'intermediaire respectivement d'un 

20 dispositlf de calcul de FFT 8i d 8j en serie avec un registre d decalage 
(9i § 9j). Dans ce mode de realisation est alors calcuie dans un mul- 
tlplicateur 10 d partir des signaux de test xj^®^ quantifie par un quantifi- 
cateur 13 couple en sortie du circuit sommateur 3. La difference 
^k'^cic-ai^ est calcuiee par un circuit soustracteur 1 1 . Un dispositlf 

25 de mise d jour 12 assure la mise d jour des gains gO) en fonction de 

Hr, Ck. 
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REVEIMDICATIONS 

1. Proc6d6 d'6galisation non lin^aire d'un signal multi-porteuse 
dans une liaison satellite caract§ris6 en ce qu'il consiste : 

5 - a effectuer (1) une post-distorsion du signal provenant du 

satellite 

- ^ decoder (2) les points de la constellation du signal nnulti- 
porteuse repu 

et § corriger (4) chaque point de la constellation par Elimination 
10 des Eohos fixes pour estimer les points de la constellation qui sent les 
plus proches des points corrig6s. 

2. Proc6d6 selon la revendication 1 caract6ris6 en ce que la post- 
distorsion du signal prend en consideration la r6ponse en amplitude du 
satellite en approximant la caract6ristique amplitude phase du signal sor- 

15 tant du satellite relativement au signal entrant dans le satellite par des 
gains complexes g^iK 

3. Proc§d6 selon la revendication 2 caract6ris6 en ce que le d6co- 
dage des points de la constellation de symboles regus a lieu en effec- 
tuant un calcul de transformSe de Fourier rapide (4 ; 8i - 8\) sur le signal 

20 post-distordu. 

4. Proc6d6 selon la revendication 3 caractdrlsE en ce qu'ii consiste 
d caiculer des coefficients de correction des raies du spectre du signal 
transform6e en effectuant une Emission pEriodique d'un paquet de signal 
test bien connu du rEcepteur, en effectuant une transform6e de Fourier 

25 de chaque paquet test regu et en la comparant (10) raie par raie d la 
transform§e de Fourier du paquet test correspondent 6mls. 

5. Dispositif pour la mise en oeuvre du proc6d6 selon Tune quel- 
conque des revendications 1^4 caract6ris6 en ce qu'il comprend un 
premier quantificateur (1) coupl6 a des amplificateurs de gains complexes 

30 (2j), un circuit sommateur (3) coupl6 sur ses entries d'opErande S la 
sortie des amplificateurs de gain complexes (2j), un dispositif de calcul 
de FFT (4) couplE entre la sortie du circuit sommateur (3) et une pre- 
miere entrEe d'op6rande d'un circuit multiplieur (6) recevant sur une 
deuxifeme entr6e d'opErande les coefficients de correction Hk §labor6 d 
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partir de remission p^riodique des paquets de test, ainsi qu'un deuxifeme 
quantificateur (7) coupl6 d la sortie du circuit multiplieur (6) pour fournir 
les points de la constellation estimSe. 

6. Dispositif selon la revendicatlon 5 caract6ris6 en ce qu'il com- 
prend en plus un dispositif de mise d jour (12) des gains complexes g(j). 

7. Dispositif pour la mise en oeuvre du proc6d6 selon Tune quel- 
conque des revendications 1^4 caract^risd en ce quMI comprend cou- 
ples entre un premier quantificateur (1) et un deuxi6me quantificateur (7) 
un circuit sommateur (3) couple par ses entries d*op6randes respectives 
d un niveau du premier quantificateur par rinterm6diaire d'un dispositif 
de calcul FFT (8-| ... Bp et d'un amplificateur de gain complexe g^lK la 
sortie du circuit sommateur 6tant reli6e au deuxi6me quantificateur (7) 
par rinterm6diaire d'un circuit multiplieur (6) recevant sur une deuxifeme 
entr6e de correction les coefficients de correction 6labor6s d partir de 
remission p6riodique de paquets de test, la sortie du deuxifeme quan- 
tificateur (7) fournissant les points a^ de la constellation estim6e. 



1/2 



2698504 



Postdistorsion 

des 
echantilLons 
de signal 



Recuperation 

des points 

de la 
constellation 



Elimination 

des 
echos fixes 



Estimation 

des 
points les 
plus proches 



FIG.1 



2698504 




FIG.3 




REPUBLIQUE FRANCAISE 



INSTITUT NATIONAL 
de la 

PROPRIETE INDUSTRIELLE 



RAPPORT DE RECHERCHE 

etabli sur la base des dernieres revendicaiions 
deposees avant le commencement de la recherche 



2698504 

N*^ d'cmcgistrenieDt 
natJODal 

FR 9214082 
FA 483852 



Catcgorie 



DOCUMENTS CONSIDERES COMME PERTINENTS 



A.D 



A,0 



Citation da document avec indications en cas de besoini 
des parties pertinentes 



Revendicaitions 
oonoeiDtes 
de la demande 



2nd European Conference on Satellite 
Communications, 22-24 October 1991, Liege, 
BE; ESA, Noordwijk, NL, 1991; 
pages 179 - 184, Rapp: "Effects of 
HPA - nonlinearity on a 4-DPSK/OFDM - 
signal for a digital sound broadcasting 
system" 

* abrege * 

* page 184, colonne de gauche, alinea 1 * 

IEEE TRANSACTIONS ON COMMUNICATIONS 
vol. 30, no. 5, Mai 1982. NEW YORK US 
pages 1233 - 1242 

SALEH ' Intermodulation analysis of FDMA 
satellite systems employing compensated 
and uncompensated TWT's* 

* abrege * 

IEEE TRANSACTIONS ON COMMUNICATIONS 
vol. 29, no. 11, Novembre 1981, NEW YORK 
US 

pages 1715 - 1720 

SALEH 'Frequency - independent and 
frequency - dependent nonlinear models of 
TWT amplifiers' 

* abrege * 

* page 1719, colonne de gauche, alinea 4 - 
colonne de droite, alinea 2 * 

US-A-4 500 984 (SHIMBO ET AL.) 

* abrege; figures 1,2 * 

* colonne 1. ligne 7 - ligne 20 * 

* colonne 1. ligne 59 - ligne 68 * 

EP-A-0 365 431 (THOMSON - CSF) 

* abrege * 

* page 7, ligne 12 - ligne 20 * 



Dale iTacUvemt de la recbeicfac 

16 AOUT 1993 



3-5,7 



DOMAINCS TECHNIQUES 
RECHERCHES QhU 0.5 ) 



H04L 
H04J 



SCRIVEN P. 



CATEGORIE DES DOCUMENTS CITES 

X : partioiUeremait pertinent a lui sail 

Y : particulierement pertinent en combinaison avec un 

autre document de la menie categorie 
A : pertinent a i'enconCre d'au moins une revendication 

ou arriire-plan technologtque gteeial 
O : divulgation non-ecrite 

P ; document Intercalaire 



T : tbforie ou priodpe i la base de finvcDtion 

E : document de brevet b^efidaot d'une date ant&ienre 

d la date de depOt et qui n'a et^ publie qn'i cette date 

dc depot ou qn'i one date postirieure 
D : cite dans la dcmande 
L : ettk pour d'auttes raisons 



&Tinembre de ia mbne famillei document conespondant 



